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I.  ACCOMPLISHMENTS 


Included  below  is  a  brief  summary  of  the  accomplishments  under  Grant 
AF0SR-80-0176. 


A.  Dual  Variable  Method  and  Navier-Stokes  Problems 

Many  natural  implicit  discretizations  of  the  Navier-Stokes  equations 
can,  with  the  proper  identification,  be  regarded  as  systems  defining  flows 
on  associated  networks.  A  set  of  network  variables  (the  "dual  variables") 
is  introduced  which  significantly  reduces  the  size  of  the  original  system 
and  hence,  economizes  on  its  solution.  The  method  avoids  the  need  to  compute 
pressures,  and  produces  velocities  that  are  exactly  discretely  divergence 
free.  Details  for  transient  problems  can  be  found  in: 


•  "An  Application  of  Network  Theory  to  the  Solution  of  Implicit 
Navier-Stokes  Difference  Equations",  R.  Amit,  C.  A.  Hall  and 
T.  A.  Porsching,  J.  Comp.  Physics,  40,  1981,  pp.  183-201. 

•  "DUVAL:  A  Computer  Program  for  the  Implicit  Treatment  of  Two- 
Dimension,  Two-Phase  Fluid  Transients",  C.  A.  Hall  and  T.  A. 
Porsching,  ICMA  Technical  Report  No.  81-25. 

•  "Numerical  Solution  of  Navier-Stokes  Problems  by  the  Dual 
Variable  Method",  C.  A.  Hall,  ICMA  Technical  Report  No.  82-42. 
Invited  Talk  at  Mathematical  Association  of  America  Meeting, 
Toronto,  August,  1982.  To  appear:  SIAD. 


The  extension  of  the  dual  variable  method  to  finite  element  discreti¬ 
zations  was  also  investigated.  Such  discretizations  of  the  Navier-Stokes 
problem  yield  systems  of  nonlinear  ordinary  differential  equations  with 
constraints  of  the  generic  form 


M  ^  +  Q  ( U )  U  +  AT  P  =  B 
AU  =  S 

for  velocities  U  and  pressures  P.  The  dual  variable  method  consists  of 
the  following  steps: 

(i)  Construct  a  basis  for  the  null  space  of  A  and 

define  the  matrix  C  =  :  •  •  •  :  $.]. 

( i i )  Construct  a  particular  solution  Up  to  (2)  then 

u  =  up  +  c  Y 

for  some  vector  y  (the  dual  variables). 

(iii )  Multiply  (1)  by  CT  (this  eliminates  P  from  (1))  and  use 
(3)  to  get 

CT  M  C  |  +  CT  Q(U)C  y  -  B'. 

(iv)  Solve  (4)  for  the  dual  variables  y  and  use  (3)  and  (1)  to 
recover  the  velocities  0  and  pressures  P. 

Research  centered  on  the  investigation  of  NS(A),  the  null  space  of  the 
finite  element  analogue  A  of  the  divergence  operator.  What  is  d,  the 
dimension  of  NS(A)?  How  can  a  basis  for  NS(A),  i.e.  C,  be  efficiently 
constructed? 

If  the  above  scheme  is  to  be  computationally  efficient  then  the  con¬ 
struction  of  C  must  be  cheap.  In  the  context  of  finite  differences,  the 


y, 

S  . 
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matrix  C  was  the  fundamental  matrix  of  an  associated  network  and  its 
"construction"  was  trivial. 

The  finite  element  investigation  was  carried  out  in  the  context  of  the 
4-node  quadrilateral  element  for  approximating  velocity  and  the  constant 
quadrilateral  element  for  approximating  pressures. 

For  virtually  arbitrary  unions  of  rectangular  elements  an  algorithm 
was  devised  for  the  construction  of  C  which  does  not  require  the  solution 
of  any  linear  system.  The  number  of  columns  of  C,  i.e.  the  dim  NS (A),  is 
determined  by  the  following  theorem: 

Theorem:  Let  A  be  the  discrete  divergence  operator  associated  with  the 
combination  of  4-node  bilinear  velocity  elements  and  constant  pressure 
elements.  If  the  domain  is  admissibly  decomposed  into  a  union  of  N 
rectangular  elements  and  there  are  M  interior  nodes  then: 

(i)  A  is  N  x  2M, 

(ii)  dim  NS(A)  =  2M  -  N  +  2. 

Note  that  (1)  and  (2)  is  a  system  of  2M  ODES  with  N  algebraic  constraints 
while  (4)  is  a  system  of  2M  -  N  +  2  (~  N)  ODES. 

For  virtually  arbitrary  unions  of  quadrilaterals  another  algorithm  has 
been  devised  for  the  construction  of  C  which  requires  the  solution  of  a 
sequence  of  relatively  small  linear  systems.  The  work  is  contained  in: 

•  "A  Finite  Element  Implementation  of  the  Dual  Variable  Method  for 
the  Navier-Stokes  Equations",  Ph.D.  Thesis  -  Frank  R.  Sledge, 

August,  1983.  ICMA  Technical  Report  83-60,  May,  1983. 
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•  "The  Dual  Variable  Method  for  Finite  Element  Discretizations  of 
Navier-Stokes  Equations",  C.  A.  Hall,  J.  S.  Peterson,  T.  A. 

Porsching  and  F.  R.  Sledge,  ICMA  Technical  Report  83-64,  October, 

1983.  To  appear:  Int.  J.  Numer.  Meth.  Eng. 

The  dual  variable  method  has  also  been  extended  to  transient  compressible 
fluid  flow  problems.  The  complication  now  is  the  introduction  of  the 
pressure  ir.to  the  continuity  equation  (density  is  now  a  function  of  pressure). 
It  is  still  possible  to  regard  implicit  finite  difference  discretizations 
of  compressible  flow  problems  as  systems  defining  "flows"  on  networks.  A 
set  of  network  variables  (dual  variables)  is  introduced  which  reduces,  by 
a  factor  of  two,  the  size  of  the  linear  system  which  must  be  solved  at  each 
time  step.  A  report  on  this  extension  authored  by  J.  Burkardt,  C.  Hall  and 
T.  Porsching  will  be  issued  in  the  near  future.  See  also 

•  "A  Numerical  Method  for  Compressible  Flows  in  Networks",  T.  A. 
Porsching,  Electric  Power  Research  Institute,  Palo  Alto,  CA,  to 
be  issued. 


Here,  as  well  as  in  many  other  applications,  F  is  a  nonlinear  Fredholm 
mapping  from  a  Banach  space  X  =  Z  *  A  to  a  Banach  space  Y,  Z  represents 
a  state  space  (usually  infinite-dimensional),  and  A  an  m-dimensional 
parameter  space.  In  the  case  of  flow  problems,  F  is  generated  by  the 
Navier-Stokes  equations  and  the  given  boundary  conditions,  z  represents 


the  state  variables  characterizing  the  flow  problem  (usually  the  components 
of  the  velocity  vector,  the  pressure,  and  so  on),  and  A  is  an  m-dimensional 
vector  of  intrinsic  parameters  (for  instance,  the  Reynolds  number  or  any 
other  dimensionless  numbers  characterizing  the  flow). 

For  practical  flow  problems,  it  is  rarely  of  interest  to  compute  only 
some  solutions  of  (5)  for  a  few  given  parameter  vectors.  Instead,  interest 
centers  on  such  concerns  as  the  variation  of  solutions  under  specified 
changes  of  the  parameters  and  their  stability  properties.  For  this  purpose, 
the  equilibrium  states  of  the  flow  problem  should  be  considered  as  points 
(z,A)  in  the  product  space  X  =  Z  x  A,  and  the  set  of  solutions  (z,A)  e  Z  x  A 
of  (5)  constitutes  an  m-dimensional  manifold  in  Z  *  A.  The  problem  then 
is  to  analyze  computationally  the  characteristic  features  of  this  equilibrium 
manifold. 

Work  on  this  topic  has  addressed  three  principal  subareas.  The  first 
involves  a  comprehensive  study  of  the  discretization  errors  arising  when 
the  original  problem  (5)  is  replaced  by  an  equation  between  finite-dimensional 
spaces.  A  theory  of  a  priori  discretization  errors  for  problems  of  the  form 
(5),  which  generalizes  the  results  of  previous  studies  in  two  important  ways, 
has  been  developed.  In  particular,  the  theory  covers  a  broader  class  of 
problems  than  the  mildly  nonlinear  problems  to  which  other  approaches  are 
restricted,  and  the  results  are  valid  for  the  general  case  of  an  m-dimensional 
parameter  space  for  any  m  >_  1 .  Details  of  this  on-going  study  appear  in 
the  following  publications: 

•  "On  the  Discretization  Error  of  Parametrized  Nonlinear  Equations", 

J.  P.  Fink  and  W.  C.  Rheinboldt,  SIAM  J.  Numer.  Anal.  20  (1983), 

732-746. 


•  "Solution  Manifolds  and  Submanifolds  of  Parametrized  Equations  and 
Their  Discretizations  Errors",  J.  P.  Fink  and  W.  C.  Rheinboldt, 

Univ.  of  Pittsburgh,  Inst,  for  Comp.  Math,  and  Applications, 

Technical  Report  ICMA-83-59,  June,  (1983),  to  appear  in  Numer.  Math. 

•  "Some  Analytic  Techniques  for  Parametrized  Nonlinear  Equations  and 
Their  Discretizations",  J.  P.  Fink  and  W.  C.  Rheinboldt,  in  Numerical 
Methods ,  Proceedings  of  the  International  Workshop  Held  at  Caracas, 
June  14-18,  1982,  Lecture  Notes  in  Mathematics,  vol .  1005,  ed.  by 

B.  Pereyra  and  A.  Reinoza,  Springer-Verlag  (1983),  108-119. 


The  second  subarea  concerns  a  study  of  computational  methods  for  analyzing 
solution  manifolds  of  problems  of  the  form  (5)  and  the  development  of 
library-grade  software  for  such  methods.  Included  in  this  study  are  the 
application  of  the  continuation  package  PITCON,  developed  by  W.  Rheinboldt 
and  J.  Burkardt  at  ICMA,  to  steady-state  flow  problems  and  the  use  of 
techniques  for  reducing  the  sizes  of  the  systems  arising  in  typical  flow 
problems.  Special  emphasis  was  placed  on  the  reduced-basis  technique. 

Details  of  this  work  may  be  found  in  the  following  publications: 


•  "On  the  Error  Behavior  of  the  Reduced-Basis  Technique  for  Nonlinear 
Finite-Element  Approximations",  J.  P.  Fink  and  W.  C.  Rheinboldt, 
Zeit.  f.  Angew.  Math,  und  Mech.  63  (1983),  21-28. 

•  "High  Reynolds-Number  Solutions  for  Incompressible  Viscous  Flow 
Using  the  Reduced-Basis  Technique",  J.  S.  Peterson,  Univ.  of 
Pittsburgh,  Inst,  for  Comp.  Math,  and  Applications,  Technical 
Report  ICMA-83-49,  (1983). 


•  "The  Reduced-Basis  Method  for  Discrete  Systems,  Part  I:  Error 
Estimates",  T.  A.  Porsching,  Univ.  of  Pittsburgh,  Inst,  for  Comp. 
Math,  and  Applications,  Technical  Report  ICMA-83-57,  (1983). 

•  "Local  Error  Estimates  for  Parametrized  Nonlinear  Equations",  J.  P. 
Fink  and  W.  C.  Rheinboldt,  Univ.  of  Pittsburgh,  Inst,  for  Comp. 
Math,  and  Applications,  Technical  Report  ICMA-84-70,  (1984), 
submitted  to  SIAM  J.  Numer.  Anal. 


The  third  subarea  consists  of  a  study  of  the  computational  detection 


and  analysis  of  singularities  on  solution  manifolds.  Such  singularities 
include  bifurcation  phenomena  of  practical  flow  problems.  Based  on  the 
general  theory  of  discretization  mentioned  above,  augmented  equations  and 
computationally  accessible  coordinate  systems  on  solution  manifolds  of 
parametrized  equations  have  been  studied  with  special  emphasis  on 
singularities.  Computational  algorithms  for  the  efficient  dete  ~n  of 
singularities  on  a  solution  manifold  are  also  being  investigat  The 
development  of  these  algorithms  stems  from  earlier  work  by  W.  .’’^boldt. 
The  results  of  this  study  are  contained  in  the  following  publications: 


"The  Role  of  the  Tangent  Map  in  Analyzing  Bifurcation  Behavior", 

J.  P.  Fink  and  W.  C.  Rheinboldt,  Univ.  of  Pittsburgh,  Inst,  for 
Comp.  Math,  and  Applications,  Technical  Report  ICMA-83-63,  (1983), 
to  appear  in  Zeit  f.  Angew.  Math,  und  Mech. 


"Folds  on  the  Solution  Manifold  of  a  Parametrized  Equation",  J.  P. 
Fink  and  W.  C.  Rheinboldt,  Univ.  of  Pittsburgh,  Inst,  for  Comp. 
Math,  and  Applications,  Technical  Report  ICMA-84-73,  (1984), 
submitted  to  SIAM  J.  Numer.  Anal. 


Included  here  for  reference  are  some  earlier  publications  on  this 


"A  Comparison  of  Methods  for  Determining  Turning  Points  of  Non¬ 
linear  Equations",  R.  G.  Melhem  and  W.  C.  Rheinboldt,  Computing  29 
(1982),  201-226. 

"Computation  of  Critical  Boundaries  on  Equilibrium  Manifolds", 

W.  C.  Rheinboldt,  SIAM  J.  Numer.  Anal .  19  (1982),  653-669. 

"Numerical  Methods  for  a  Class  of  Finite-Dimensional  Bifurcation 
Problems",  W.  C.  Rheinboldt,  SIAM  J.  Numer.  Anal .  15  (1978),  1-11. 

"A  Locally  Parametrized  Continuation  Process",  W.  C.  Rheinboldt  and 
J.  V.  Burkardt,  ACM  Trans,  on  Math.  Software  9  (1983),  215-235. 
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•  "Algorithm  596:  A  Program  for  a  Locally  Parametrized  Continuation 
Process",  W.  C.  Rheinboldt  and  J.  V.  Burkardt,  ACM  Trans,  on  Math. 
Software  9  (1983),  236-241. 

C.  Mixed  Finite  Element  Methods  for  Stationary  Navi er- Stokes  Problems 

Conforming  mixed  finite  element  methods  for  incompressible  viscous  flow 
problems  were  examined.  Preliminary  work  consisted  of  studying  linearized 
Navier-Stokes  equations  such  as  the  Oseen  equations.  This  work  culminated 
in  a  report  entitled 

•  "On  Conforming  Mixed  Finite  Element  Methods  for  Incompressible 
Viscous  Flow  Problems",  M.  D.  Gunzburger  (University  of  Tennessee), 

R.  A.  Nicolaides  (Carnegie-Mel 1  on  University)  and  J.  S.  Peterson 
(University  of  Pittsburgh),  Comp.  &  Maths,  with  Appls.  8  (1982), 

167. 

In  this  paper  the  asymptotic  rates  of  convergence  are  given  for  four  finite 
element  discretizations. 

This  work  has  been  extended  to  the  inhomogeneous  stationary  (nonlinear) 
Navier-Stokes  equation  with  inhomogeneous  boundary  data.  Specifically,  the 
problem  considered  is  to  find  u  and  p  satisfying 

9u 

-vAu  +  u.  t—  +  grad  p  =  f  in  n  (6) 

J  dx. 

div  £  =  g  in  ft  (7) 

U|r  =  g_  (8) 

for  given  f  e  H~^(ft),  g  e  L2(ft)  and  g_  e.  H^(fi)  such  that 
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f  g  =  f  g.-n 
Jn  J  r 

where  r  is  the  boundary  of  ft.  Here  Hr(ft)  for  r  >  0  denotes  the 

r-th  order  Sobolev  space  associated  with  the  set  ft;  H”r(ft)  denotes  the 

dual  space  of  H  (ft);  and  [^(ft)  represents  the  space  of  functions  which 

are  square  integrable  in  ft.  For  the  particular  weak  formulation  of 

(6)-(8)  considered,  it  is  possible  to  show  existence  and  uniqueness  of  the 

solution  under  certain  conditions  on  the  data.  In  addition,  optimal  error 

estimates  for  the  velocity  measured  in  the  (ft)  norm  and  for  the  pressure 
2 

measured  in  the  L  (ft)  norm  for  the  Galerkin  approximation  to  the  weak 
problem  are  derived.  Since  the  discretization  process  leads  to  a  system  of 
nonlinear  algebraic  equations,  iterative  techniques  (simple  iteration,  Newton 
iteration  and  Modi fi ed- Newton  iteration)  are  used  to  solve  this  system. 

Rates  of  convergence  for  each  scheme  are  derived.  Numerical  examples  which 
use  an  efficient  choice  of  finite  element  spaces  are  also  provided.  This 
work  is  presented  in  the  following  report: 

•  "On  Conforming  Mixed  Finite  Element  Methods  for  the  Inhomogeneous 
Stationary  Navier-Stokes  Equations",  J.  S.  Peterson  and  M.  Gunzburger 
(Carnegie-Mellon  University),  Numer.  Math.  42  (1983),  173-194. 


An  additional  area  of  interest  for  the  Navier-Stokes  problem  is  to 

2 

determine  the  error  in  the  velocity  measured  in  the  L  (ft)-norm.  To  obtain 

such  a  result  it  was  first  necessary  to  improve  upon  the  existing  theory  of 

the  approximation  of  boundary  value  problems  with  inhomogeneous  essential 

boundary  conditions  such  as  the  inhomogeneous  Dirichlet  problem  for  second 

order  elliptic  partial  differential  equations.  An  optimal  error  obtained 
o 

in  the  L  (ft)-norm  was  then  obtained  for  the  (nonlinear)  stationary  Navier- 
Stokes  equations.  This  analysis  of  the  finite  element  method  for  problems 
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with  inhomogeneous  essential  boundary  conditions,  along  with  numerical 
examples  for  problems  posed  on  polyhedral  domains,  is  presented  in  the 
report 

•  "On  Finite  Element  Approximation  of  Problems  Having  Inhomogeneous 
Essential  Boundary  Conditions",  J.  S.  Peterson  and  G.  Fix 
(Carnegi e-Mel  Ion  University)  and  M.  Gunzburger  (Carnegie-Mel Ion 
University),  Comp.  &  Maths,  with  Appls.  9  (1983),  687. 

The  work  described  above  deals  with  the  Navier-Stokes  equations  in 
primitive  variable  formulation.  Also  of  interest  is  the  streamfunction- 
vorticity  formulation  of  the  Stokes  and  Navier-Stokes  equations.  Finite 
element  algorithms,  have  been  derived  which  require  low  continuity  finite 
element  spaces  and  do  not  require  any  artificial  specification  of  the 
vorticity  at  solid  boundaries.  In  addition,  methods  for  handling  multiply 
connected  domains  have  been  analyzed  both  theoretically  and  computationally. 
This  work  is  presented  in  the  paper 

•  "On  the  Finite  Element  Approximation  of  the  Streamfunction-Vorticity 
Equations",  M.  Gunzburger  (Carnegie-Mel Ion  University)  and  J.  Peterson, 
Proc.  5th  IMACS  International  Symposium  on  Computer  Methods  for  Partial 
Differential  Equations,  1984. 

D.  Krzhivitski -Ladyzhenskaya  Finite  Difference  Equations 

In  1966,  Krzhivitski  and  Ladyzhenskaya  presented  a  finite  difference 
discretization  of  the  Navier-Stokes  problem  for  which  they  prove  unconditional 
stability  in  the  discrete  L2-norm.  Their  scheme  also  generates  approxi¬ 
mations  which  converge  to  a  weak  solution  of  the  continuous  problem  as  the 
discretization  parameters  tend  to  zero.  Their  analysis  required  homogeneous 
boundary  data. 

Under  this  grant  the  analysis  of  the  K  &  L  difference  scheme  was 


extended  in  two  respects: 


(i)  Existence  and  uniqueness  theorems  for  solutions  of  the  stationary 
K  &  L  equations  were  established  for  both  homogeneous  and  inhomogeneous 
boundary  conditions.  The  results  for  the  homogeneous  boundary  conditions 
were  announced  in  the  paper  and  Ph.D.  thesis 


"On  the  Krzhivitski-Ladyzhenskaya  Difference  Method  for  the 
Stationary  Navier-Stokes  Equations",  A.  Cha-Lin  and  T.  A. 
Porsching,  Proc.  10th  IMACS  World  Congress,  Montreal,  Canada, 
Vol .  2,  R.  Peskin,  ed.,  (1982). 
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•  "On  the  Krzhivitski-Ladyzhenskaya  Difference  Equations  for  the 
Stationary  Navier-Stokes  Equations",  Ph.D.  Thesis  -  Anna  Cha-Lin, 
December,  (1982),  Univ.  of  Pittsburgh,  Inst,  for  Comp.  Math,  and 
Applications,  Technical  Report  ICMA-82-50,  January,  (1983). 

(ii)  As  a  step  towards  handling  curved  domains  consideration  was  given 
to  extensions  of  the  K  &  L  difference  scheme  to  a  differential  system  that 
arises  when  the  Navier-Stokes  system  is  transformed  through  the  introduction 
of  a  curvilinear  coordinate  system.  This  required  the  derivation  and 
analysis  of  a  new  class  of  finite  difference  approximations.  A  new  7  point 
stencil  was  derived  and  was  proven  to  be  unconditionally  stable.  For  a  given 
mesh  gauge  h,  a  piecewise  bilinear  function  u^  was  constructed  which 
interpolates  to  the  finite  difference  solution,  and  it  has  been  shown  that 
as  h  -*■  0,  u^  converges  to  the  weak  solution  of  the  continuous  problem. 

Details  of  this  work  are  contained  in: 


•  "Analysis  of  Finite  Difference  Discretizations  of  Navier-Stokes 
Problems  on  Mapped  Domains",  Ph.D.  Thesis  -  John  Ellison,  August, 
(1982),  Univ.  of  Pittsburgh,  Inst,  for  Comp.  Math,  and  Appli¬ 
cations,  Technical  Report  ICMA-82-41 ,  July,  (1982). 


E.  Binary  Gas  Mixture  Flow  Through  Combustors 


In  an  attempt  to  reduce  the  development  cycle  costs  associated  with 
design  of  gas  turbine  engine  combustors,  mathematical  combustor  models  are 
being  employed  to  provide  information  about  performance  trends  and  to 
predict  velocity,  pressure  and  thermodynamic  property  profiles  in  simulated 
practical  combustion  environments.  It  has  been  demonstrated  that  the  dual 
variable  method  can  be  applied  to  the  predictive  model  of  the  fluid  dynamics 
associated  with  an  axially  symmetric  centerbody  combustor  being  studied  at 
WPAFB.  This  was  discussed  in  the  report 

•  "Non- Isothermal  Flow  Through  an  Axially  Symmetric  Centerbody 
Combustor  via  the  Dual  Variable  Method",  C.  A.  Hall  and  T.  A. 
Porsching,  Univ.  of  Pittsburgh,  Inst,  for  Comp.  Math,  and  Appli¬ 
cations,  Technical  Report  ICMA-82-38,  May,  (1982). 

A  cooperative  effort  between  ICMA  and  the  Aero  Propulsion  Laboratory, 
WPAFB  is  underway  to  incorporate  the  dual  variable  reduction  technique  into 
their  combustion  modelling  program  TEACH-T. 


F.  Numerical  Solution  of  Convection- Conduction  Problems 

The  convection-conduction  equation  arises  in  fluid  dynamics  as  the 
vorticity  transport  equation.  In  convection  dominated  situations  (high 
Reynolds  number)  numerical  solutions  may  contain  severe  nonphysical 
oscillations.  Upwind  finite  difference  methods  and  finite  element  methods 
have  been  developed  over  the  years  to  prevent  such  spatial  oscillations. 

Upwind  finite  element  schemes  were  developed  under  this  grant  which 
lead  to  diagonally  dominant  systems  of  positive  type: 


v. 


•  "Numerical  Solution  of  Convection-Conduction  Problems  by  Nonconforming 
Upwinded  Finite  Elements",  I.  Christie,  Univ.  of  Pittsburgh,  Inst,  for 
Comp.  Math,  and  Applications,  Technical  Report  ICMA-82-43,  August, 
(1982). 


This  guarantees,  for  example,  that  the  solution  satisfies  a  maximum  principle, 
a  property  which  is  not  shared  by  all  finite  element  discretizations  as  shown 
in 


•  "The  Maximum  Principle  for  Bilinear  clements",  I.  Christie  and  C.  Hall 
Univ.  of  Pittsburgh,  Inst,  for  Comp.  Math,  and  Applications,  Technical 
Report  ICMA-82-46,  August,  (1982) ,  I.  J.  Num.  Meth.  Enq.  20,  (1984), 
549-553. 

•  "Maximum  Principles  and  Finite  Elements",  I.  Christie,  Proc.  Ninth 

S.  African  Symposium  on  Numerical  Mathematics,  J.  Snyman,  ed.,  Durban, 
July  18-20,  (1983). 
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